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ABSTRACT 

Cataclysmic variables undergo periodic nova explosions during which a finite mass 
of material is expelled on a short timescale. The system widens and, as a result, the 
mass-transfer rate drops. This state of hibernation may account for the variety of 
cataclysmic variable types observed in systems of similar mass and period. In the light 
of recent changes to the theory of nova ignition and magnetic braking we investigate 
whether hibernation remains a viable mechanism for creating cataclysmic variable 
diversity. We model the ratio of time spent as dwarf novae (DNe) to nova-like systems 
(NLs). Above a critical mass-transfer rate the system is NL and below it a DN. The 
dominant loss of angular momentum is by magnetic braking but the rate is uncertain. 
It is also uncertain what fraction of the mass accreted is expelled during the novae. 
We compare the models of the ratios against the period of the system for different 
magnetic braking rates and different ejected masses with the ratio of the number of 
observed NLs to DNe. We deduce that a rate of angular momentum loss a factor of 
ten smaller than that traditionally assumed is necessary if hibernation is to account 
for the observed ratios. 
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1 INTRODUCTION 



Cataclysmic variables JWarnerl fl9951 are interacting close 
binary stars in which a white dwarf is accreting from a com- 
panion. In the vast majority of cases the companion is a low- 
mass main-sequence star or red dwarf and the periods are 
measured in hours. The red dwarf is transferring hydrogen- 
rich material to the white dwarf by Roche lobe overflow at 
a rate, 10 -11 to 10 -8 MQyr -1 . As the hydrogen settles on 
the surface, the base of the layer becomes degenerate and 
heats. When enough, 1CP 5 to 1O~ 3 M0, has accumulated it 
ignites. The degenerate conditions lead to a thermonuclear 
runaway and a nova explosion. The system brightens by ten 
or so magnitudes and most of the layer if not all of it to- 
gether with some of the underlying white dwarf material is 
blown off. 

Though some of this material may fall back on to the 
white dwarf and some of it may be accreted back by the red 
dwarf, a significant fraction escapes carrying with it only the 
intrinsic angular momentum of the white dwarf. The system 
expands such that aM = const, where a is the separation 
and M the total mass which has just decreased. If the red 
dwarf had no atmosphere, that is a distinct boundary just 
filling its Roche lobe, it would no longer do so and mass 
transfer would cease. This is the idea of hibernation first 
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documented bv lShara et alJ l)l986l) . At that time they were 
motivated by the apparent lack of cataclysmic variables in 
the solar neighbourhood when compared with theoretical 
estimates and the frequency of classical novae in M31. The 
complete detachment during hibernation would render the 
majority of cataclysmic variables unobservable at any given 
time. 

Because the red dwarf has a slightly extended atmo- 
sphere it actually must overfill its Roche lobe in order to 
maintain an equilibrium mass-transfer rate. The ratio of the 
atmosphere scale height to the radius of the star is compa- 
rable to the ratio of the mass lost from the system to its 
total mass, so the change in separation is not sufficient to 
shut off mass transfer altogether but does reduce the rate 
significantly. This change in rate could be responsible for 
the diversity of cataclysmic variable types at a given period: 
dwarf novae that undergo disc instability outbursts require 
a lower mean mass-transfer r ate than nova-likes that do 
not s uffer the disc instability dFaulkner. Lin fc Papaloizoul 
1983). However the idea lost favour when it was realised 
that the nova ejecta would temporarily form a common en- 
velope around the system and consequential angular mo- 
mentum l oss would actually cause low er-period systems 
to shrin k ifL^io^iorar ie. fc Ritterl 11993) . Recent work by 
lNelemarrTet^lTl200(lH ndicated that the common-envelope 
phase is very poorly understood and that in some cases sys- 
tems with tenuous envelopes can actually expand. Since we 
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wish to compare the simplest case we neglect any common- 
envelope phase in this work. 

A detailed a nalysis of the period distribution was made 
bv lShafterl (I1992T) . He emphasized that the lack of dwarf no- 
vae between 3 and 4 hr is difficult to reproduce theoretically 
with any reasonable magnetic braking law. He speculated 
that the solution might lie in a better understanding of any 
of additional angular momentum loss alongside mass trans- 
fer, weak but unnoticed magnetic fields that disrupt the in- 
ner parts of the disc, correlation between white dwarf mass 
and orbital period or hibernation. Among these he deduced 
that hibernation might prove the most promising and it is 
this that we explore here in more detail. 

We investigate the hibernation scenario for cataclysmic 
variable diversity in the light of two recent changes to the 
accepted model of cataclysmic variables. These are 

(i) the relation between the mass of hydrogen that must 
be accreted to ignite Mj gn and the secular accretion rate 
jTownslev fc Bildstenll2004h 

(ii) the possibility that magnetic braking, that drives the 
evolution of cataclysmic v ariables, could be much weaker 
than previously thought iAndronov. Pinsonneault fc Sills! 
2003). 



We consider hibernation in its simplest form to deduce what 
combinations of J, Mi gn and M e j, the mass actually ejected 
in the nova, could be consistent with the observed variety of 
cataclysmic variable types. 



2 SIMPLE ANALYTIC MODEL 

We first consider a simple analytical model to get some idea 
of what fraction of the time is spent in hibernation and 
then consider the simplest numerical model that includes 
the pertinent effects of the atmosphere and variations of 



2.1 Binary Systems 

The cataclysmic binaries we consider consist of a white 
dwarf, of mass Mi, and a low-mass main-sequence star, of 
mass M 2 and radius R 2 . 

The hydrostatic and thermal equilibrium radius of the 
main-sequence star can be approximated by, 



Rq 



M 2 



(1) 



The Roche lobe radius 7?li is the radius of a sphere 
which encloses the same volume as that enclosed by the 
last stable equipotential surface around star i. The Roche 
lobe radius of star 2, for a restricted range of mass ratio 
q = Ma/Mi, is 



Rl2 _ _2_ /MA 1/3 
a ~ 3 4 / 3 \M ) ' 



< q < 0.8, 



where M = M i + M 2 is the total mass of the system 
JPaczvnskilll97ll) . The simplicity of this formula makes it 
useful for analytic work within the restricted but useful 
range of q. When R 2 > Rl 2 mass is transferred from the 
red dwarf to the white dwarf. 

Kepler's third law, that relates the angular velocity Q 



of each star and the binary system to the total mass M and 
the separation of the two stars a, is 

2 GM 



(3) 



When the red dwarf is filling its Roche lobe, R 2 = Rl2, we 
combine equations (1), (2) and (3) to find 

3/2 

/' - — — I — ^ I Ma, (4) 



Rq 

(2G)V2 

where P = 2ir/Q is the orbital period. So in this model 
the orbital period depends only on the mass of the main- 
sequence star when R 2 = Rl 2 - 

We assume the spin angular momentum of the stars 
is negligible compared to the orbital angular momentum. 
Hence the total angular momentum of the system is 

J = Mia?H + M 2 a\VL, (5) 

where a\ and a 2 are the distances from star 1 and star 2 
respectively to the centre of mass of the system. We can 
rearrange 

aiMi = a 2 M 2 and a = ax + a 2 (6) 
to find 



M 2 



(7) 



and similarly for a 2 . Substituting into equation (5) we have 

(8) 



MiM 2 a 
J — — — — a il 



M 

Using Kepler's law, equation (3), and equation (8) we can 
rearrange to get the separation 

_ J 2 M 

a ~ GMIMI 1 > 

and the angular velocity 
G 2 MfMi 



n = 



J 3 M 



(10) 



in terms of the three independent variables, J, Mi and M 2 



2.2 Nova Eruption 

When a critical mass Am, which we examine later, has accu- 
mulated on the surface of the white dwarf star, thermonu- 
clear reactions ignite in the degenerate material. Suppose 
the entire layer of Am is expelled from the system in a nova 
explosion. Since AMi = AM — — Am, the angular momen- 
tum of the system changes by 



A J = AMiajn = - Am 



/aM 2 



£1 



(11) 



(2) 9± = 



We want to find an expression for the change in separation 
during the nova eruption. Equations (8) and (11) give 

A J _ AmM 2 

T ~ ~ 



(12) 



MMi 

Differentiating equation (8) gives 

AJ _ A Mi AM _ AM Aa AO. 

~T-^iT + ^4 2 W + a + 1T' {S) 

but, because M 2 is constant during the eruption, AA1 2 = 0. 
Differentiating equation (3) we get, 
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(14) 



(15) 



AO, _ AM Aa 
~Q~ ~ ~W ~ ~o~" 

We can combine equations (12), (13) and (14) to get 
the familiar result 
Aa Am 
~a~ = ~W' 

Because of the mass loss from the system, the separation 
of the stars increases during the explosion. Equations (12) 
and (15) ensure that a(Mi + M 2 ) = const during the erup- 
tion. The Roche lobe radius must also increase as the separa- 
tion increases. Differentiating equation (2) and substituting 
equation (15) we find 

A7? L2 4Am 

-R^ = W (16) 

After the eruption R2 < Rl2 and so the mass-transfer rate 
decreases sharply, switches off in this model. In this low 
mass-transfer rate the cataclysmic binary is said to be in a 
state of hibernation. 



2.3 Conservative Case 

We now consider a fully conservative case for which M 
const and J = so 



J = Mi M 2 a 

J Mi M 2 2a ~ ' 

Since Mi = —M2 we can substitute to find 

-=2(?-l)# 
a yH 'Mi 

and 

P 
P 



3(9- 



M2 
M2 



(17) 



(18) 



(19) 



where q = M2/M1. These two equations show that, as the 
mass of the red dwarf decreases, a and P only decrease if 
q > 1. But because the Roche lobe continues to shrink as q 
decreases, mass transfer doesn't end at q = 1. Differentiating 
equation (3) we get 



RL2 
i?L2 

and combining with equation (18) we find 



a M 2 
a + 3M? 



= Ml ( 9 - !\ 
fl L2 M 2 \ q 3J 



(20) 



(21) 



This shows that the minimum Roche lobe radius is when 
q — 5/6. But for most observed cataclysmic binaries q < 5/6 
and steady mass transfer occurs. For a red dwarf to be in 
steady contact with its Roche lobe there must be something 
causing loss of orbital angular momentum. 

So after an eruption the stars gradually get closer to- 
gether because of the loss of angular momentum which 
brings the main-sequence star back into contact with its 
Roche lobe. The mass-transfer rate increases once more and 
the cycle is repeated. We have many eruptions. 

We want to find the range of mass ratios q for which 
mass transfer is stable. Since we have, by differentiating 
equation (1) 



R2 
R2 



M2 
M 2 



(22) 



and in order to be stable we need 

Rl2 R2 
Rl2 > R2 ' 



(23) 



by comparing equations (21), (22) and (23) we therefore 
must have 



q < 4/3. 



2.4 Analytic Ratio of Times 



(24) 



After a nova eruption the main-sequence star no longer fills 
its Roche lobe because the orbit has expanded. Here we 
assume that the rate of angular momentum loss |J| remains 
constant until the next nova explosion. 

When the system is detached (when there is no mass 
transfer) Mi = AI2 = M = 0. So from one nova explosion, 
the time until mass transfer begins again is the time it takes 
until Rl2 shrinks back to R 2 - From equations (13), (14) 
and (16), 



AJ 
J 



\J\t d 



Aa 
~2a~ 



AR L 



2Am 



J 2a 2R L 2 3M 
So the time the system spends detached is, 
J 2Am 

U = 



j\ 3M 



(25) 



(26) 



When the system is semi-detached mass is transferred 
from the red dwarf to the white dwarf, Mi = —M% and 
M = 0. The Roche lobe of the red dwarf is filled so that 
R2 — Rl2 and R2 — Rl2 so 
Ai?L2 Aa AM 2 AR 2 AM 2 



_ Aa AM 2 
i?L2 ~ ~a~ + SM 2 ~ R2 ~ ~M~2 
Rearranging this we find 

Aa _ 2AM 2 
~~a ~ 3M 2 

Substituting (28) into (13) we get 
\J\t s Am Am Am 



AJ 
J 



Am Am 



J Mi M 2 3M 2 

The time spent semi-detached is 
3q) 



J . (4 
-Am- — — — 

\J\ 3M 2 



(27) 



(28) 



(29) 



(30) 



This model predicts the ratio of time spent detached to the 
time spent semi-detached to be, 



(31) 



U = (l + g)(4-3g) 
U 2q 

We plot it in Fig. Q for later comparison with our numer- 
ic al model and note th at this is consistent with the result 
of lLivio fc Sharal (119871) that hibernation should be deepest 
when q is close to one. 



3 OBSERVED CATACLYSMIC VARIABLES 

Classical novae are a class of novae and cataclysmic vari- 
ables that have had only a single observed eruption. Their 
brightness ranges from 6 to greater than 19 magnitudes. 
Dwarf novae (DNe) are a class of cataclysmic variables 
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Figure 1. The solid line is t s /id against the mass ratio q (equa- 
tion 31). We compare with the dashed line which is the ratio of 
the time spent as a nova-like to dwarf nova against q with a con- 
stant period P = 3hr, J factor of 0.1 and M; gn 
numerical model (section 5). 



where the prime distinguishes this from a prescription that 
includes the effect of varying the mass-transfer rate. 

However the white d warf is heated by accretion s o the 
rate also depends on /M^. lTownslev fc Bildstenl <l2004l) have 
calculated the critical mass which causes ignition on the 
white dwarf. We make a linear fit to their figure 8 and find 



(M) 



M ign _ 

M ~" ^lO^Moyr- 1 
where 

-4 /M< 



a = 3.2 x 10" 
and 



Ui J 



1.231 



= 0.534 



Mi 
Mo 



(34) 



(35) 



(36) 



When Mi gn has accumulated on the white dwarf we have a 
nova explosion. We consider three different models, two in 
which the mass ejected is constant in all the explosions, 

(i) M e j = const = 1 x 1O" 3 M 

(ii) M ej = const = 5 x 1O~ 4 M 

and one in which the mass ejected in the nova equals the 
amount of mass needed for ignition, 

(iii) Mej = M ign . 



that have had several observed outbursts. These outbursts 
however range in brightness from 2 to 5 magnitudes. They 
are associated with a disc instability. Nova-like systems 
(NLs) include all non-outbursting cataclysmic variables. 
Their spectra indicate that they could be novae in a pre- 
or post-eruption stage. 

The nature of the cataclysmic variable, DN or NL, 
should depend only on the mass-transfer rate which itself 
should only depend on the system masses and period. How- 
ever at any given period some of each type are seen. 



4 NUMERICAL MODEL 

We now estimate the radius and l uminosity o f the m ain- 
sequence star with the formulae of iTout et al. (1996) and 
use th e more accurate formula for the Roche lobe radius for 
star i teggletonlll983Tl . 



0.49g,; 



2/3 



0.6 g f 3 +log(l + ^ /3 )' 



(32) 



where q; is the mass ratio Mi/Ms-i. This formula is valid 
for all q. 



4.1 Ignition Mass 

By considering cold white dwarf models iFuiimotol dl982l) 
calculated the critical pressure at th e base of the hydrogen 



rich layer for ignition. Based on this lTruran fc Livid iff?. 
derived 



M[ gn _ 9.42 x 10- 
M Q _ Mi/M G 



(i?i/lO 9 i? )' 



(33) 



4.2 Rate of loss of Mass 

When the main-sequence star is overfilling its Roche 
lobe mass is transferred on to the white dwarf so that 



Mi = -M 2 



(37) 



We approximate the atmosphere of the red dwarf by a simple 
isothermal model. The rate of mass loss from the red dwarf 
is then given by, 



M 2 



-Ms 



(38) 



where AR is the amount by which the star overfills its Roche 
lobe, 



AR = R 2 - i? L2 . 



(39) 



When the red dwarf underfills its Roche lobe, R2 < Rl2, AR 
is negative and the mass-transfer rate is very small. When 
the star overfills its Roche lobe, R2 > Rl2, AR > and the 
mass transfer rate increases. 

The scale height H of the unperturbed isothermal at- 
mosphere of the main-sequence star is 

p(r) 



H = 



dp(r) 
dr 



(40) 



so we need to find p(r), the density at position r. The equa- 
tions of an isothermal atmosphere are 

^W=-p(r) fl and P(r) = ^p(r), (41) 
dr p 

where P(r) is the pressure, 5R is the gas constant, T c is the 
effective temperature at the surface of the main-sequence 
star and p is the mean molecular weight. The atmosphere is 
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a mixture of hydrogen, helium and metals in various states 
of ionization and in molecules. So we use p, equal to 1. The 
gravitational field strength at the surface of the red dwarf 
is, 



9 = 



GM 2 
R\ 



(42) 



The surface temperature of the main-sequence star is found 
by Stefan's law, 

1/4 



T e = 



( L 2 



Integrating the equations in I41H for p, 

p = Po e 



(43) 



(44) 



where p is some constant. Hence we find the scale height 



H 



(45) 



The factor in front of the exponential in equation (38) is an 
estimated maximum mass-transfer rate equal to the mass of 
the secondary star divided by its dynamical timescale. We 
have verified that varying it by large factors has almost no 
effect on our results because AR simply adjusts to compen- 
sate when mass transfer is stable. 



4.3 Rate of Loss of Angular Momentum 

Angular momentum is continually lost from the system. For 
the closest systems gravitational radiation is the most im- 
portant cause of mass transfer. The rate of loss of angular 
momentum in gravit ational radiation from two point masses 
in a circular orbit is jLandau fc Lifshitzlll95lh 



Jgr _ 32G 3 M 1 M 2 (M 1 +M 2 ) 
J ~ 5c 5 a 4 



(46) 



In wider systems however magnetic braking is the more 
important mechanism. The rate of change of angular mo- 
mentum ow ing to magnetic b r aking is uncertain. Based on 
the wo rk of ISkumanichl il972l) . iRappaport. Verbunt fc Joss! 
l|l983l) postulated a rate of the form 

(47) 



Jmb = -5.83 x lO~ 16 (fl 2 /-R0) 3 (r2yr) 3 M Q 7?|yr" 



However recent work by 

lAndronov. Pinsonneault fc Sills! <l2003l) based on the 
spins of stars in clusters suggests that this rate is orders 
of magnitude too high. But such a low rate would not be 
able to drive the interrupted magnetic braking model for 
the period gap between two and three hours in cataclysmic 
variables. In this model a system evolves to lower period as 
M 2 drops at such a high rate that the red dwarf is out of 
thermal equilibrium. At a period of three hours, a mass of 
about O.3A40, magnetic braking abruptly falls off, perhaps 
because the red dwarf has become fully convective. The 
mass-transfer rate drops and the red dwarf shrinks back 
to its equilibrium radius and must wait for some much 
weaker angular momentum loss mechanism to shrink the 
orbit until Rl 2 is again reduced to R 2 . With the high rate 
the period gap is very well reproduced. With the low rate 
the stars are not far enough from thermal equilibrium and, 
even if magnetic braking were interrupted, the gap would 



be much too narrow. If such a low rate proves to be correct 
an alternative reason for the period gap must be found. 

4.4 Irradiation 



The nova decay time is 50 to 100 yr dLivio fc Sharall987t) . So 
although irradiation of the red dwarf from the white dwarf 
may increase the mass-transfer rate initially, it will only be 
for a small fraction of the inter nova period of several thou- 
sand years. Though we find increased mass transfer while the 
white dwarf is hot the time for which it is is negligible and 
we ignore it. We note th at the irradiation cycles discussed by 
iBiining fc RitteJ i2004l) operate on a much longer timescale 
and would have to survive several episodes of hibernation. 



5 RATIO OF NUMBER OF OBSERVED NLS 
TO DNE 

We select the s ystems from Ritter's Catalogue 
jRitter fc Kolbl 120031) which are classified as nova-likes 
or as dwarf novae and divide them into period bins with at 
least five systems of each kind in each bin. We then divide 
these numbers to find the observed ratio of NLs to DNe. 
Though there are obvious selection effects, such as the fact 
that dwarf novae, owing to their variability, are far more 
easily spotted, we can hope to reproduce trends in this 
ratio. 

Fig. |21 shows that the ratio peaks at a period of just 
over 3hr then decays as the period increases further. After 
a period of 5 hr there are very few observed binaries but the 
data do suggest that the ratio increases again slightly with 
increasing period. We use the period of the peak ratio 3.5 hr 
and the fact that it has dropped by a factor of twelve by 
4.4 hr to test our models. We do not expect to reproduce 
the exact height of the peak because different selection ef- 
fects apply to DNe and NLs. We exclude identified polars 
and intermediate polars but we cannot be certain that oth- 
ers, particularly IPs, do not contaminate the sample of NLs. 
However their contribution is small, certainly when com- 
pared with other probable selection effects, at the relatively 
high periods we consider here. 

5.1 Critical Mass- Transfer Rate 

We calculate the times spent as a dwarf nova toN and as a 
nova-like £nl by considering the mass-transfer rate. Fig. |3] 
shows an example of the mass-loss rate of the red dwarf for 
three nova explosions. 

There is a critical mass-transfer rate above and below 
which the system appears to be in two different states. This 
critical rate of mass loss from the main-sequence star is 



M 2 



-8.08 xlO 1 



( aH 
V0.3 



(l + .) 7/8 ( 



p. 7/4 

hr/ 



- X (48) 



iFaulkner. Lin fc Papaloizoul ll983). where an = 0.1. Dwarf 
novae lie below the critical value. They have (M 2 ) < M 2cT i t 
whereas nova-likes all lie above the critical M 2cl a relation- 
ship. 

Just after a nova explosion, the mass-loss rate from the 
red dwarf is very small. The time from the explosion until 
the mass-loss rate has increased up to the critical rate is 
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Figure 2. The ratio of the number of observed nova-likes to dwarf 
novae against the period. Each point is a ratio based on at least 
five of each type of system. Poisson noise error bars are therefore 
typically a factor of two of the height of the point. 



Figure 3. The mass-loss rate against time for a system with 
initial mass ratio q = 0.8 and period P = 3.5 hr with M e j = Mi gn 
and J = O.Umb m equation 1471 . The dashed line is the critical 
rate below which the system would undergo dwarf nova outbursts 
and above which it would not. 



tDN. The binary then appears to be a nova-like from this 
time until the next eruption, tnh- We calculate the ratio of 
these times. 

Ratio = — . (49) 

If the critical mass-transfer rate is higher than the rate 
of the model for the whole cycle, then we are predicting that 
there are no nova- likes, only dwarf novae. So we have tNL = 
toN = time between nova eruptions. Similarly, if the critical 
rate is lower than the mass-transfer rate, we have only nova- 
likes, tDN = and £nl = time between nova eruptions. 

In Fig. we compare this ratio for P — 3 hr against 
q with the simple analytic calculation in section 2.4. At a 
period of 3 hr there are both NLs and DNe at approximately 
all mass ratios and we see that both models are indeed very 
similar. 



Table 1. Models with different factor in j and M e \ 



J Factor M e j/M0 Peak Peak Ratio after peak 

Period/hr Ratio atP = 4.4hr 



1.0 M ign 4.5-7.25 inf 

1.0 5 x 10" 4 6.0 inf 

1.0 1 x 10" 3 3.5 4.7 4.3 

0.5 M ign 5.75 13.5 

0.2 M ign 3.25 3.1 2.0 

0.15 M ign 3.25 2.8 1.0 

0.1 Mign 3.0 2.4 0.2 

0.1 5 x 10~ 4 <2.0 >4.0 0.2 

0.1 1 x 10~ 3 <2.0 >3.4 0.2 

0.1 M! <3 >1.2 0.6 



5.2 Model Ratio Averaging 

In order to compare our models with the observed data we 
need to plot the predicted ratio against the period but the 
ratio also depends on q. For a given period, we average the 
ratio over varying q. Averaging by adding the total times, 
£nl and tuN, over a flat q distribution and then finding the 
ratio of these total times is our preferred method. It auto- 
matically takes account of system lifetimes. The longer the 
time, the more weight it is given in the averaging. In the 
observed data, the longer a system is in a state, the more 
likely it is to be observed. We do not try to model the birth 
rate of various systems with population synthesis because 
we consider this to be too uncertain. We therefore implicitly 
assume all cataclysmic variables equally likely to form. 



6 RESULTS 

To examine the effects of a much lower angular momen- 
tum loss rate we multiply the formula for magnetic braking 
(equation 47) by a constant factor. We try combinations of 
different J factors and M e j and compare the ratios in Ta- 
ble 1. 

In Table 1 we list various models and properties which 
can be compared with the observations. As we have dis- 
cussed, selection effects make any truly quantitative com- 
parison difficult but we have identified key features in Fig. 
that ought to be reproduced. We also concentrate on peri- 
ods above the gap because we are interested in what can be 
said about magnetic braking irrespective of whether or not 
it is interrupted. Then we see a peak in the ratio of NLs to 
DNe just above 3 hr with a rapid decrease to higher periods. 
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P/hr 

Figure 4. The ratio against period for the models with J factor = 
0.1,0.15,0.2 (from bottom to top) and M e j = M- lgn . 

By 4.4 hr it should have fallen by a factor of 12. A record of 
'inf for the peak ratio indicates no dwarf nova e. 

With the full lRappaport. Verbunt fc Jossl (Il985t brak- 
ing rate (J factor = 1) we cannot reproduce these features. 
When just the accreted mass is blown off there are no DNe at 
all between about 4.5 and 7.25 hr. If a constant 5 x 10 _4 Mq 
mass is ejected then there is a peak in the ratio at 6 hr. Only 
when M j = 1O _3 M0 do we find a peak at about 3.5 hr but 
then the ratio has hardly fallen b y 4.4 hr. In any ca se such 
a large ejected mass is unrealistic. iKolb et al.l i200lf) inves- 
tigated the effects of hibernation at these standard rates of 
magnetic braking. They deduced that M G j > 2 x 10 -4 Mq 
would be necessary to account for the observed spread in M 
at a given period above 3 hr. 

It is not until we have reduced the braking rate by a 
factor of ten that we can reproduce the observed features. 
There is little variation between the three different M e j , all of 
which show a peak in the ratio at about 3 hr with a drop by 
a factor of five or more by 4.4 hr. We also made one model 
with the old cold white dwarf ignition mass M c j = M[ gn . 
Here, the ratio is only 1.2 at P — 3hr predicting about 
equal numbers of NLs and DNe. Because DNe should be 
easier to observe than NLs this is inconsistent with the data 
too. 

Fig.^Jshows the ratios as a function of period for J fac- 
tors of 0.1, 0.15 and 0.2. All have the correct shape but as the 
J factor increases the curves flatten, no longer reproducing 
the drop of a factor of 12 by 4.4 hr. 

For a J factor of 0.01, which gives a rate close to that 
of lAndronov. Pinsonneault fc Silla (l2003f) . all of the mass- 
transfer rates are lower than the critical value and so we 
see only dwarf novae. The only way this could be consis- 
tent with observations would be if another mechanism such 
as irradiation or magnetic cycles can periodically raise the 
mass transfer by an order of magnitude. This would then be 
responsible for the variety of cataclysmic variable types. 



7 CONCLUSIONS 

For the full lRappaport. Verbunt fc Jossl (I1983T) rate ( j fac- 
tor of 1) only the unrealistically high M c j of 10 _3 Mq follows 
the observed trend but then the variation of ratio with pe- 
riod is too small. 

With J reduced by a factor of ten all the models show 
the right trends with the ratio rising to a peak at a pe- 
riod of about 3hr and then falling distinctly by a period 
of 4.4 hr. We cannot then distinguish between the differ- 
ent mass ejection models but do find a distinct improve- 
ment when M e j = Mi gn when pro per account is taken 
of wh ite dwarf heating by accretion jTownslev fc Bildstenl 
120041) . Hibernation is most effective at reproducing observed 
trends when our J factor is a bout 0.1. This is ten times 
smalle r than that proposed bv lRappaport. Verbunt fc Jossl 
jl983j) and ten times larger than that proposed by 
lAndronov. Pinsonneault fc Sills! J2003ft . 
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